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Abstract 

The lowest order pomeron loop is calculated for the leading conformal weight with full 
dependence of the triple pomeron vertex on intermediate conformal weights. The loop is 
found to be convergent. Its contribution to the pomeron Green function begins to dominate 
already at rapidities 10-^15. The pomeron pole renormalization is found to be quite small 
due to a rapid fall of the triple pomeron vertex with rising conformal weights. 

1 Introduction 

In the framework of QCD, in the limit of large number of colours, strong hadronic interactions 
are mediated by the exchange of BFKL pomerons which split and fuse by triple pomeron 
vertices. This picture can be conveniently described by an effective nonlocal quantum field 
theory [1]. A remarkable property of this theory is its inherent conformal invariance, which 
is broken by interactions with colliding hadrons. In terms of Feynman diagrams contributions 
standardly separate into tree diagrams and diagrams with pomeron loops. For reactions with 
heavy nuclei the tree diagram contribution is enhanced by factor A 1 /^ for each interaction 
and so dominates. This dominating part can be summed by the Balitski-Kovchegov equation 
for DIS on a heavy nucleus [2, 3] or by a pair of equations constructed by the author for 
nucleus- nucleus scattering [4]. However reactions with single hadrons require taking into 
account also loop diagrams. There has been several attempts to calculate the contribution 
of a single pomeron loop [5, 6] with a crude approximation for the triple pomeron vertex 
and contradicting results. In [6] it was found that the magnitude of the loop is so small 
that it gives no significant contribution up to extaordinary high energies (rapidities of the 
order of 40). Lately there were several claims to sum all loop contributions in the colour 
dipole approach or in the so-called reaction-diffusion formulation of the scattering mechanism 
[7, 8, 9, 10, 11]. However very crude approximations made from the start do not allow to 
consider these results even minimally reliable. These circumstances give us a motivation to 
reconsider the contribution of a single pomeron loop by the conformal invariant technique 
developed in [1]. This technique in fact greatly simplifies the derivation and allows to fix 
numerical coefficients, uncertainties in which in our opinion were one of the reason why the 
results of [5] and [6] turned out to be different. Most important we also use the exact form 
of the triple pomeron vertex, which appears very different from its approximate value used 
in the previous calculations. 

Our results first demonstrate that the pomeron loop, with all contributions taken into 
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account, is finite and does not need renormalization, in contrast to the old local Regge-Gribov 
model. Its numerical value is found to be small indeed, but not so small as calculated in [6]. 
As a result, its influence becomes visible at much lower energies than claimed there. With 
realistic values for the QCD (fixed) coupling constant its contribution starts dominating 
already at rapidities y ~ 10 — 15. This of course means that taking loops into account for 
reactions with single hadrons is necessary already at present energies. 

The paper is organized as follows. In the next section we recall some elements of the 
conformal technique introduced in [1] to be used in loop calculations. Using this technique 
we calculate the loop contribution in the next section. Section 4. presents our numerical 
results and discusses influence of the loop contribution on the pomeron Green function. 
Some conclusions are drawn in the last section. Technical details and comparison with [6] 
are discussed in three appendices. 



2 The pomeron interaction diagrams 

Feynman diagrams for the pomeron interaction are built from the pomeron propagator and 
triple pomeron vertex. The pomeron propagator 

where 1 = {r^\r^} are the initial coordinates of the two reggeized gluons, 2 = {r^\r^} 
are their final coordinates and y\ — y 2 is the rapidity difference, satisfies the equation 

(^ + H)g y _ y ,(l\2))=5(y-y')V- l 2 V 2 H{l\2), (2) 

where H is the BFKL Hamiltonian [12]. The triple pomeron vertex can be read off the 
interaction Lagrangian 

Li = / I 2 \ - 1 )<P 2 )L 12 tf 3 + h.c, 3 

where <j> and $ are the two fields which describe the propagating pomerons, operatorLi2 is 
defined as 

L 12 = r\ 2 V\vl (4) 

and the fields in (3) are to be taken at the same rapidity. 

In absence of interaction with external hadrons the theory is conformal invariant. This 
makes it convenient to pass to the conformal basis formed by functions (in complex notation) 
[12] 

Here n = {n, v, ro} = {h, ro}, h = (1 — n)/2 + iu, h = (1 + n)/2 + iu, with n integer, v real 
and two-dimensional transverse ro, enumerate the basis. In the following, for clarity, we shall 
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sometimes write h as a set of two numbers {n, v}. We also pass from rapidity y to complex 
angular momentum j = 1 + u: 

9y = / . (6) 

J a—too fc 

Then the propagator can be presented as 

g u (l\2) = '£E li (l)E;(2)g u , h , (7) 

M>0 



where 



/U n=— co 

with 



°° /-CO 1 /• 

E = E / <*V / d 2 ro, (8) 
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a h = a^ = — 9 ■ Z3 - T - A . (9) 



2 z^ 2 + n 2 /4' 
The conformal propagator is 

11 

9u>,h = -j , (10) 

L nu UJ — LOh 

where uih are the BFKL levels 

u h = 2a(-0(l) -ReV>(fr))> ® = ^ £ (H) 

and 

4tt 8 

k = • (12) 

Q , n-\-\,vQ , n—\,v 

The triple pomeron vertex can be presented in the conformal basis as 

r(l|2,3)= E ^i(l)^ 2 (2)^ 3 ( 3 ) r Mi|M2M3' (13) 

Ml>A t 2A 1 3>0 

The dependence on the intermediate cm. coordinates R±, R2 and R3 is fixed by conformal 
invariance 

r Mi|M2M3 = ^12 12 - R 23 3 - R 31 1 - ( a - /■) - r fti|fc 2 ^3- ( 14 ) 

Here (a.f.) means the antiholomorhic factor. Powers are known 

11/ XV 

«12 = -- + - \ni - ni - n 3 ) + i(y x - v 2 + 1/3), 

11/ X ■/ 

«23 = + 2 ( ni + 71,2 + ~ H^ 1 + ^2 + ^3), 

a 3 i = -- + -(n 3 - ni - n 2 ) + + z/ 2 - ^3) (15) 

The powers in the antiholomorhic factor a's are obtained by changing rii — > — nj. In the 
lowest approximation the conformal vertex rj^^ ^ 3 = h 2 h 3 w ^h h = h* was introduced 
and studied by Korchemsky [13]. 
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Using (7), (13) and the orthonormlization property of states (5) one can perform inte- 
grations over gluon coordinates and substitute them by summations over conformal weights 
and integration over center-of-mass coordinates, as indicated in (8). With the known expres- 
sions for the propagator and vertex given by (10) and (14) respectively, one can then write 
expressions for any Feynman diagram directly in the conformal basis. One should only take 
into account that 'energies' to are conserved at each vertex. 



3 The pomeron self-mass 

3.1 The pomeron full Green function 

The full pomeron Green function G w (l|2) can also be written in the form similar to (7): 

G„{l\2) = Y J E»(l)E;{2)G^ h . (16) 

The Schwinger-Dyson equation expresses the full Green function via the pomeron self-mass: 

Gu > h = TTn TWy — • ( 17 ) 

where X^/j is the pomeron self-mass in the conformal basis. 

Similar to (7) the pomeron self-mass in the gluon coordinate space can be written as a 
sum over conformal eigenstates 

E(l|l')= E ^^(1)^(1'), ( 18 ) 

Mi ,Mi 

where the self-mass in the conformal representation is 

_ 8^Nl f dui ^ (o) ( 

MlMi ^2 / 2-Ki ^ Mi|M2,M3 LT ^ LT ^ i M2M3|M'i- V ly J 
M2.M3 

and the suppressed dependence on to is obvious from its conservation at the vertexes. From 
its conformal invariance it follows that 

= <W'i S Mi> ( 20 ) 

where 

<W = Snn'${v - v')$ 2 {roo')ah- (21) 

Summation over ^ an d H3 in (19) includes integration over two center-of-mass coordinates 
i?2 and R3. Dependence on them comes only from the vertex functions and indicated in (14) 
and (15). 

The second vertex part T^^^i is defined by the expansion 

r(2,3|l)= £ ^w(i)^(2)£ w (3)r 

M2M3|m'i - 

(22) 

Mi>M2,M3>0 



We use the property of functions 

J d(l)E*, (1)^(1) =5^, »,fi'>0, 
where d(l) = d 2 rid 2 r[/rf v , to find 

WaK = / d(l)d(2)d(3)^ 1 (l)^ 2 (2)^ 3 (3)r(2,3|l). 
However G(2,3|l) = G(l|2,3) and 

E li l ) = E A l )i P- = v(n^ -n,v ^ -v), 

so that 

WaK = / d(l)d(2)d(3)E* Pl (l)^ 2 (2)^ 3 (3)r(l|2,3) = 1^,^. 
So we find an integral in (19) 

— i rfijj j2p r?ai2 e>c*23 r? a 3l 7?* "12 r>* "23 r?* "31 p/323 r>* o* $23 o* 

— / « -K3-K12 ^23 ^31 ^12 -^23 ^31 ^1'2 ^23 ^31' ^1'2 ^23 ^31' 

The full set of powers is 

a 12 = -- + -(n 2 - ni - n 3 ) + - z^ 2 + 1/3), 

«23 = + ^(ni + n 2 + ra 3 ) + - v 2 - vz), 

q 3 i = -- + -(n 3 - ni - n 2 ) + + f 2 - ^3) 
Q12 = + ^( _n2 + n i + n 3) + ^1 - ^2 + ^3), 

"23 = ~ + 2"( _ni ~ n 2 - "-3) - + ^2 + ^3), 

«31 = - ^ + \(~ n ^ + n l + n 2) + +"2- V 3 ) 

Pl'2 = ~2 + ^( _n2 + n 'l + n3 ) + ^"^l + U2 ~ ^3). 

#23 = + ^(~ n 'i ~ n 2 - n 3 ) + + z/ 2 + ^3), 
1 1 

031' = - 2 + 2 ( _n3 + n 'l + + - ^2 + ^3) 

/3i'2 = + 2^2 _ n 'i ~~ n 3) + + vz), 

023 = + ^( n 'l + n 2 + n 3) + + ^2 + ^3), 
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for = ~\ + \( n z - n'i - n 2 ) + i{-v[ - "2 + i*) (28) 



So we find 

S mK = 8C ^2 Vc J r h?|/ Hl ,hs G «i,ft2 G w-«i 1 fts r feh3|V 1 ^i.^.^i.^ / i)- ( 29 ) 

According to conformal invariance property (20) this expression has to be proportional to 
Spin' and this has to be valid with any values for Gh and r^j^g, since with any values for 
these quantities conformal invariance of the Green and vertex functions is fulfilled. Taking 
unity for the vertex functions and deltas for the Green functions in the conformal basis we 
find that the integral (27) itself has to be proportional to S^^i and independent of R\ 

I(R 1 ,X 1 MX) = 8 ni ^5fa-rt)P(Ru>)F(h 1 \h2M- (30) 
The pomeron self-mass is expressed via F according to (20) and (21): 

h,2,h,3 

3.2 The lowest order pomeron self mass in the conformal representation 

Calculation of F(h\\h 2 hs) is described in Appendix 1. It is found that 

F(h 1 \h 2 ,h 3 ) = a hl . (32) 

This gives for the self-mass 

ggfog f dui sr r(°) 



OftjSh! - ^ 2 c J -^t ^tl^hfi^iMG^-^i^Shih^Fih^h^). (31) 



_ Sa^NZ f duj! ^ ( o) , . 

^wM~ ^2 J 2-ni ^ hi\h2,h 3 ^uiM ^u-u-lM L h 1 \h 2 ,h 3 - \ M ) 
Here 

n=-oc> u \ / 

In the lowest order we substitute the Green functions G by pomeron propagators g given 
by (10) and the full vertex T by its lowest order expression T^°\ Introducing the explicit 
expression for 1^ and doing the integral over u)\ we then obtain 



ni,Ti2=— oo 

where 



a A N 2 +oo oo \ T hLh 2 \ 2 

S in c Y] / dvidv 2 — b hl b h2 , (35) 



h = K ' v = [!/2 + 1)2/4]^ + („_ ^2/4]- ( 36 ) 

One observes that the contributions from conformal weights n\ = ±1 and n 2 = ±1 
seemingly diverge at small z^i and v 2 , when in the denominators we find factors v\ or v\. 
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However one can demonstrate that at least for symmetric state {n, v} (that is with even n) 
the vertex cannot be coupled to any of antisymmetric states {ni, u±} and {712^2} (that is with 
odd n\ or 722) (see Appendix 2.), so that this divergence is in fact absent. The situation for 
antisymmetric initial state {n, v} is not clear, since the 3-pomeron vertex was derived only 
for a symmetric state (gg-loop). In the following we assume the initial state {n, v\ to be 
symmetric (n even). 

In fact the lowest order vertex function in the conformal representation h ^ = Q(h, ^1,^2) 
was studied in [13]. It was found to be highly complicated. It was expressed in [13] in terms 
of the Meier function G44. In [13, 14] only its value for the leading conformal weights 
h = hi = h-2 = 1/2 was found: 

r ( °\ , = \, \) = 7766.679. (37) 
2 ' 2 ' 2 v z z z 

Self-mass (35) is obviously an analytic function of uj with a left-hand cut. For a particular 
term in the sum over n\ and n 2 the cut goes from uj ni! o + Wti 2 ,o to —00. The rightmost cut 
corresponds to ri\ = n-2 = and starts at uj = 2A where A = 4a In 2 is the BFKL intercept. 
The discontinuity of £ across the cut is given by 



a 4jy2 too roo (Q) 

= ~ i ^^~ ^ J dVldV2 \ T h\h±M\ K^i^ ~ -Wfca). (38) 
ni,n2=— oo ^ 

It is negative imaginary and different from zero for uj < 2aj]y 2 = 8(5 In 2. Therefore the original 
leading BFKL pole at uj = uJh=i/2 acquires an imaginary part and splits into two complex 
conjugate poles which remain on the physical sheet of the complex u-plane in contrast to 
normal theories where the pole goes under the cut onto the unphysical sheet. 

It is not difficult to find the asymptotic behaviour of the Green function in the conformal 
representation as a function of rapidity 

Gy,h = J ^ y " G u,h- (39) 



At y — > 00 it is dominated by the contribution from the rightmost cut extending from 
uj = 2uji/ 2 = 2 A to uj = —00 

G y h [° dwe^Disc 1 , - =— [° duje^- DiSCS ,^ , 9 . (40) 

2mJ-oo l h uj - u> h + l h ^oj,h 2mJ-oo \uj - uj h + lh^u>,h\ 

The leading contribution comes from conformal weights n\ = 712 = and small v\ and zv 2 . 
So we get 



G 



's' c 



/ duje yu) f dv\di>2v\v\ — 0u \ 0ui '°' V2 _ ! _J — ; (41) 
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where we used that at small v l nu c± 1. Expanding at small v±, u 2 in the standard manner 

^ 1(2) = A - au 2 1{2) , a = 14C(3)a (42) 
and doing the integral over uj we get 



G 



a i N c 2yA f°° j„ j ,,2,,2^-?/a(^ 2 +^) 1^0^10^1,0^2 1 

87r A 7o K 2 + V4) (^ 2 + l/4) 

Taking all factors which have finite values at v\ = v 2 = out of the integral we find finally 

The asymptotic is negative and of the order l/N 2 as expected. So in principle it belongs to 
a higher order in the expansion in 1/N C . 
For N c = 3 and taking a s = 0.2 we have 

A = 0.530, a = 3.21, (7 = 6.26. (45) 

With these values the Green function with a loop becomes greater than the bare one already 
at y ~ 10. To compare, in [6] the contribution from the pomeron self-mass was found to 
be very small due to partly a smaller numerical factor (see Appendix 3.), but mainly due 
to a different manner of studying the asymptotic. The authors of [6] assumed that the 
propagators around the loop are also in the asymptotical regime and accordingly restricted 
integration over v to small values. This introduced additional damping of the contribution 
due to weight v 2 . With their value of the coupling constant a s = 0.3 they found that the 
loop influence becomes significant only at rapidities greater than 40. However in fact the 
propagators around the loop enter not in their asymptotical regime, which greatly enhances 
the magnitude of the loop contribution. 

To see the influence of the loop on the position of the pole in the pomeron Green function 
(its 'mass renormalization') one has to calculate the self mass in the vicinity of the pole. 
Restricting to the Green function for the leading conformal weight h = 1/2 and taking into 
account in the sum (35) only the leading intermediate conformal weights with n\ = n 2 = 
one has to evaluate S Wjl ^ 2 as a function of energy uj far from the tip of the cut at oj = 2A. 
Then one has to know S7(l/2, hi, h 2 ) as a function of hi = 1/2 + iv\ and h 2 = 1/2 + iv 2 at 
v\ and v 2 greater than zero. In previous calculations [6] a very crude estimate of ^ u ,i/2 was 
obtained assuming that O, does not substantially change in the essential integration region and 
can be approximated by the known $1(1/2, 1/2, 1/2). If one uses our numerical coefficient in 
the expression for the self-mass then with this approximation the real and imaginary parts of 
^w,i/2 at u; = A turn out to be of the same order as the BFKL intercept A. The two complex 
conjugate poles corresponding to the "physical" pomeron are then found to be significantly 
different from the bare pomeron pole : 

uj$ = 0.473 ± 0.027 i (46) 
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v 1 

Figure 1: Absolute values of S7(l/2, 1/2 + iv±, 1/2 + ^2) as a function of z^i at different z^- 
Curves from top to bottom correspond to U2 = 0.0, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5. 

(recall that the "bare" intercept was real and equal to A = 0.530). The real part of the 
intercept is diminished in accordance with the conclusions in [6]) although this change is 
not too small due to a greater coefficient in our self-mass. The important fact is that the 
pole acquires an imaginary part of the same order as the change in the real part. This fact 
challenges our standard renormalization methods, since it cannot be compensated by adding 
extra terms to the original Lagrangian. 

However this calculation obviously overestimates £ at values of 00 far from uo = 2 A, since 
in fact the three pomeron vertex function O rapidly diminishes in this region. Calculation 
of 0(1/2, hi, I12) for different h\i = 1/2 + 2^1,2 requires a complicated numerical procedure. 
Wee used the formulas derived by Korchemsky in [13] in the form of integrals over variable 
x in the interval [0,1]. 1 The found 0(1/2, hi, /12) rapidly diminishes with v\ and v-i- In 
Fig.l we illustrate this behaviour showing [0(1/2, 1/2 + ivi, 1/2 + iu2)\ as a function of ^2 at 
different t>\. 

Values of £^,1/2 given by Eq. (35) with 0(1/2, hi, h2) depending on hi and hi are shown 
in Fig. 3 as a function of to. One observes that at to = A the self-mass becomes quite small. 
With a s = 0.2 

£a,i/2 = -0-0058 - i0.0569, (47) 

The renormalization of the pomeron intercept is thus insignifiant: the two complex conjugate 

1 We have checked that all of them are indeed expressed by the Meijer function as indicated in [13] except 
Ji. For the latter the expression in terms of the Meijer function probably contains a misprint. 
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Figure 2: Real and imaginary parts of S w l / 2 as a function of lo/A in the region A < lo < 2A 
for a s = 0.2, A = 0.530, with the 3-Pomeron vertex depending on internal conformal weights 

poles corresponding to the "physical" pomeron are now 

u$ = 0.524 ±0.057 i (48) 
to be compared with the estimate (46), which neglects the vertex dependence on v. 

4 Conclusions 

We have calculated the lowest order loop contribution to the pomeron Green function in the 
conformal technique. The main novelty of our calculation is the use of the triple pomeron 
vertex with full dependence on the intermediate conformal weights. On the technical side, 
our loop is found to carry a much greater numerical coefficient as compared with the old 
calculations in [6]. As a result, the loop is not at all so inocuous as stated in [6]: its 
contribution begins to dominate already at rapidities of the order 10-^15. On the other hand, 
due to the variation of the triple pomeron vertex, its influence on the pomeron pole (mass 
renormalization) is found to be quite small in agreement with [6]. 

On the theoretical side we have seen that the pomeron self mass is finite, at least in 
the lowest order, so that mass renormalization is not obligatory, unlike the old local Regge- 
Gribov pomeron model. We have also confirmed that due to the wrong sign of in front of the 
self-mass the bare pomeron pole splits into two complex conjugate ones, which stay on the 
physical sheet, contrary to what happens in the 'normal' theory. 

Finally we stress that we have limited ourselves to the lowest order loop. With a small 
coupling constant this enables us to study the asymptotic of the behaviour of the Green 
function only up to a certain finite (large) rapidity determined by the condition a s exp Ay ~ 1, 
that is 

1 . 1 
y < — In— . 

A a s 
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At larger values of y loops of higher order step in. The true asymptotic at y — > oo obviously 
requires summation of all loops. In our opinion, the path to achieve this goal is still quite 
long. 

5 Appendix 1. Calculation of F(hi\h,2, h^) 

Using representation (30) we are going to obtain F by putting n[ = n\ and integrating over 
R\. Since integral (27) depends only on the difference Rw, we can put R[ = 0. After that, 
with n[ = ri\ the integral takes the form 

I(R 1 ,n 1 ,u 1 ,u[) 

= J d 2 R 2 d 2 RzR^R%{ 1 Rl 2 &12 Rl l &zl R%i i+P2 ' i ^ ■ ( 49 ) 

Integration over R\ leads to an integral 

J(R 23 ) = J d 2 R 1 R^ 2 R 3 *f 1 R* 12 &12 R* 1 &[il = f d 2 R 1 R° 31 R^ 2 Rf 31 R* 10 ai2 , R = R 23 . (50) 

Its dependence on R 23 can be esily found by rescaling R\ = Rqz: 

J(R 23 ) = R^ 3 m+2iui Rh m+2ivi j(hi,h 2 ,h 3 ), (51) 

where 

j(h 1 \h 2 ,h 3 )= J d 2 zz a31 (l-z) ai2 z* &31 (l-z*f 12 . (52) 
Putting this result into (49) we find 

J d 2 R 1 I(R 1 ,n 1 ,v 1 ,v[) 
= j(hi\h 2 ,h 3 ) J d 2 R 2 d 2 R 3 R^fR* 23 ^ 3 R^ 2 R^'R* 2 ^ 2 Rf 3 ^, (53) 

where 

723 = -1 - ni + + vi), 723 = -1 + ni + i{v[ + v{). (54) 
To calculate (53) we use a Fourier transform: 

R a R* & = g(a, a) J ^ q -^ q *~^ e ^W)/\ ( 55 ) 

where 

g(a, a) = 2^-"2 1+ "+" ^ + ^ . (56) 

r(-a) 

The inverse transformation is 

qP q *P = gtf, 0) J d 2 RR- 1 -PR*- 1 -Pe-W R+t > R *V 2 , (57) 
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where 

~g(0, 0) = g~\-l - 0,-1- (3) = i-H)^^^^. (58) 

Of course 

f ^ e i^ R+ ^/ 2 = 5\R). (59) 

Using (55) we present products R^fR^ 23 , R^' 2 Rf^ and R% 31 ' Rf^' as Fourier trans- 
forms to obtain 

J d 2 R 1 I(R 1 ,n 1 ,v 1 ,v[) = d(h 1 ,h 2 ,h 3 ) J j[-^d 2 R 2 d 2 R^ 1 -^qr 1 - 123 

q^^ 21 ' g2 _1_/321 ' g^ 1 "^ 31 ' ql' 1 ' 1331 ' e ~ i (<ll R h+<lt R 23)/2+i(q2R* 2 +qiR2)/2-i(q3R* 3 +q* 3 R3)/2 ^ (gO) 

where 

d(hi, h 2 , h 3 ) = j{hi\h2, h 3 )g(j 23 j 23 )g(0 21 > 21 >)g(0 31 > 31 >). (61) 
Integrations over R 2 and R 3 give 

(27r) 4 ,5 2 (g 1 -( ?2 )5 2 ( (7l -g3) 

and we find 

J d 2 R 1 I(R 1 ,n 1 ,u 1 ,u[) = d(h 1 ,h 2 ,h 3 ) J ^q~^ 3 -Psv q * -fa, -fa, (62) 
We have 

-3 - 723 - 0iv ~ 031' = -3 - 723 - fhv ~ 031' = -1 - i{v\ - v[). 

The integral gives \8[y — u'), so that we are left with the calculation of factor d(h±, h 2 , h 3 ) 
defined by (61) which reduces to the calculation of j(hi\h 2 ,h 3 ) defined by (52). 

Presenting i?" 31 i?J Q31 and R^Q-R\ Q a21 in integral (50) as Fourier transforms we get 



J(Ro) = 9(a2ia 2 i)g(a 3 ia 31 ) 

Integration over R\ gives (2ir) 2 5 2 (qi — q 2 ) and we are left with 

J(Ro) = g(a 2 ia 2 i)g(a 3 ia 3 i)I 3 , (64) 

where 

h = J J^rf2 fll g 1 -2-5 2 i-fi3i (z .-2-a 21 -a3i e i( 9l fl S + g iflo)/2_ (65 ) 

According to (55)and (58) integration over q\ gives 

J 3 = R-^+ 2i ^R*^+ 2i ^g(-l - m - 2iu u -1 + m - 2iv ± ), (66) 
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which means 

j(hi\h 2 ,h 3 ) = g(a 2 ia 2 i)g(a 3 ia 3 i)g(-f23^h): 

where we used (54) . 

Collecting our results we find 

F(hi\h 2 ,h 3 ) = ^ 9(021, a 2 i)g{a 31 a 31 )g(j 23 , 723)5(723723)5(^21' hv)g{fcvhv), 
which after trivial calculations simplifies to 

F{ hl \h 2 M) = 2(^174) • 

that is to (32). 

6 Appendix 2. Coupling of the 3-pomeron vertex to pomerons 
with different parity 

In our paper [15] we considered coupling of the 3-pomeron vertex to pomeron states sym- 
metric under the interchange of gluon coordinates (of positive parity). However in loops the 
intermediate pomerons in principle may be of different parity, both positive and negative. So 
we have to know how they couple to the 3-pomeron vertex. 

To this end we have to recall that originally the vertex (coupled to the loop) consists of 
4 terms which differ by permutation of gluon momenta or coordinates, so that the vertex as 
a whole is symmetric in all 4 gluons which couple to outgoing pomerons. In the momentum 
space we had the vertex coupled to the loop as 

V(l, 2, 3, 4) = (l/2) 5 2 (G(l, 2 + 4, 3) + G(l, 2 + 3,4)+ G(2, 1 + 4, 3) + G(2, 1 + 3, 4)) . (70) 

where G was a known function. 

We couple it to two outgoing pomerons with wave functions ^i(l, 2)\l/ 2 (4, 3). In [15] 
passing to the coordinate space we in fact considered only the first term arguing that the rest 
will give the same for symmetric \t(l, 2) and ^(3, 4). Now we repeat this derivation with all 
4 terms in (70) assuming that 

*i(2,l)=Pi*(l,2), * 2 (4,3) = P 2 * 2 (3,4), P 1>2 = ±1. (71) 

We recall that the first term in (70) leads to the following triple pomeron contribution 

T (1) = ~-M f Y dy I ^fl 2 ^ *i(n, r 2 ; Y - y)* 2 (r 2 , r 3 ; Y - y)r\^\*{r^ y). 
4 4vr J Jo J r( 3 r( 2 r^ 3 

(72) 

As compared to [15] we have a factor 1/4 since there the vertex summed all 4 terms in (70) 
which gave identical contribution. To relate this formula directly to the first term in (70) 
we introduce coordinates of the 3d gluon and also of the initial pomeron into it. We also 



(67) 

(68) 
(69) 
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suppress all rapidity dependence together with integration over y which is of no importance 
for our purpose. Denoting 



44^3' r y 1 ^' ) - Jl Jl Jl 



we then have 

4 



T« = c f nd 2 r J dV 1 d 2 ^F(l,2,3)<5 2 ( 24)5 2 (l'l)<5 2 (2'3)^i(l,2)^ 2 (4,3)L i;2 ^ 3 (l / 2')- (73) 
J i=i 

In this form it is clear what we get from the rest terms in (70). 

To find the contribution from the second term in (70) we have to interchange 3^ 4 in the 
vertex (not touching the wave functions). We get 

r 4 

T(2) =c n d2r ^M^M^(l ) 2 , 4)(5 2 ( 23)(5 2 (l / l) ( 5 2 (2 / 4)^i(l,2)^ 2 (4,3)L 1 ,2^3(l'2 / ) 
i=i 

= c J d 2 r 1 d 2 r 2 d 2 r 4 F(l,2,4)* 1 (l,2)^ 2 (4,2)L 14 *3(14), (74) 
or changing integration variable r 4 — > r 3 

=c J d 2 rid 2 r 2 d 2 r 3 F(l,2,3)^i(l,2)* 2 (3,2)L 14 ^3(13) = P 2 T«. (75) 

In a similar manner we find the contribution of the third term in (70) by interchanging in 
TW l<-»2: 

T(3) =c / n d2r ^M^M^(2,4,3)<5 2 (14)«5 2 (l'2) ( 5 2 (2 / 3)^ 1 (l,2)^ 2 (4,3)L 1 , 2 ^ 3 (l'2 / ) 

i=l 

= c /" d 2 r 2 d 2 r 3a ! 2 r4F(2,4,3)^i(4,2)^ 2 (4,3)L 23 * 3 (23). (76) 
Changing integration variables 2 — > 1 and 4 — > 2 we find 

= c J d 2 nd 2 r 2 d 2 r 3 F(l,2,3)^i(2,l)^ 2 (2,3)L 23 ^ 3 (13) = PiTW. (77) 

Finally to find the last term coming from (70) we have to interchange both 1 <-> 2 and 3^4 
in T(l): 

r 4 

T^=c n^ d M^M^(l,2,4)(5 2 (13)<5 2 (l / 2)5 2 (2 / 4)vI/ 1 (l,2)^ 2 (4,3)L 1 , 2 ^ 3 (l'2 / ) 
J i=i 

= d 2 ri d 2 r 2 d 2 r 4 F{l, 2, 4)* 2 (4, 1)L 24 * 3 (24). (78) 
Changing 2 — > 1, 4 — ► 3, 1 — ► 2 we find 

T( 4 ) = c | d 2 r 1 d 2 r 2 d 2 r 3 F(l,2,3)^ 1 (2,l)* 2 (3,2)L 13 ^3(13) = PiP 2 T«. (79) 
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Summing all terms we find that the total triple pomeron contribution is 

4 

T = Y / T (i) = T^(l + P 1 + P 2 + P 1 P 2 ) =rW(l+Pi)(l + P 2 ). (80) 
i=i 

The factor multiplying TW is zero if any of outgoing pomerons is antisymmetric in its gluon 
coordinates and equal to 4 if they both are symmetric. So we find that the pomeron can 
split only into two symmetric pomerons. In this derivation it was implicitly assumed that the 
original pomeron is symmetric in its gluons (corresponding to its coupling to the qq loop). 
It is an open question to see what happens if the original pomeron is antisymmetric in its 
gluons. 

7 Appendix 3. Comparing normalization with [6](BRV) 

We first compare normalization of their 3-pomeron vertex V . From Eq. (BRV.74) we have 

V-<&M*1vAn.V, (81) 
y/N?-l 32 8 y V ' 

where V is just the Bartels vertex K 2 ^ without factors. On the other hand our definition 
starts with [15] 

ZD 2 = 2g 2 G = -2g 4 N c K 2 ^ 3 ®D 2 , (82) 
which implies that our vertex 

r = -g A N c V (83) 

(taking into account that we have to take 1/2 in view of the fact that D 2 is twice the pomeron). 
This means that our 3-pomeron vertex is related to BRV as 

vrS/2 

V = — r. (84) 

Next we compare normalizations of the Green functions and impact factors. From (BRV. 6) 
we conclude 

G BRV = 2ttG. (85) 

The amplitude is in the lowest order 

A = — [ pL-QBRV^RV 1 = !£ f f^^ 2 - ' 



2 J (2vr) 3 1 2 k 2 {q-kf 2 J (2vr) 2 1 'k 2 {q-kf 
wherefrom the relation between the impact factors is 

§ BRV = V2^$. (86) 

The BRV impact factor in the conformal representation 

= J ^E H (k, q -k)* B ™(k,q-k) = 
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= w?l a ^ e '^ t ^ I dp * m ** m (*>?m »)*(*i. fa). 

We conclude from this 

= 7^f<^ (87) 



Now we are in a position to analyze (BRV.(61)). Separating is/2 we have in the u> 
representation: 

AA= ( c 4V y 

4(2vr) 4 V16(2tt) 5 / 

J dy±dy 2 j dvdv 1 dv 2 v 1 v{vl^l m * 9Y- yu h9y 1 - V2 M9yi-V2M9y2,h = 

2ir (2vr) 2 / C 4 y \2 



-V 

(27r) 10 4(27r) 4 Vl6(2vr) 5 / 

J dyidy 2 j dvdv x dv 2 V 1 v\vl<$> h <$>h * 9Y -y 1 ,h9y 1 -y2,hi9yi-y2,h29y2,h 
where h = 1/2 + iu, /i 1 ( 2 ) = 1/2 + 2^1(2) and 

C 4V = 2 1 V/ 2 a^ Q 1 r 
Since v 2 = 7r 4 /(2a/ l ), we find 

(27r) 10 4(27r) 4 Vl6(27r) 5 / 
f f dv du x dv 2 „ 

J dy\dy 2 J '^^~^-®h®h 9Y-y 1 ,h9y 1 -y2M9yi-y2M9y2,h, 

which implies 

T BRV _ r f dl/ l d »2 „ n /oo^ 

^h, yi - y2 — ° J 2a^2a^~ yi ~ V2 ' hl yi ~ V2,h2 

where 

_ 2tt 2vr) 2 / C 4V x2_ 12 _ a^iV 2 (I 1 1\ 

(2vr) 1 4(27r) 4 Vl6(27r)5^ " vr 2 v 2'2'2y- 1 j 

Comparing with our expression we find 



*l RV = (90) 
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